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INFRARED  BdSSIVm  OF  DIAfTOKIC  CAS^  FOR 
THE  ANHABNONIC  VISITING  ROTATOR  W)IEL 

W.  Mallaaus 
A.  Thouon 

ABSTRACT 

In  order  to  conpute  the  emlsslvlty  of  a  dlatoodc  gas,  a 
siaqpllfled  model  of  a  diatomic  molecule  Is  assikmed:  an  enhar¬ 
monic  oscillator  vlth  the  first  approximation  1[o  the  vibration- 
rotation  Interaction*  For  a  given  band,  the  fj(*equ«icy  of  emitted 
radiation  Is  expressed  as  a  quadratic  function  of  the  quantum 
nuniber  m,  vhlch  Is  solved  to  express  m  as  a  fuijtctlon  of  u*  This 
expression  Is  substituted  for  m  In  the  equatloi^s  used  for  com¬ 
puting  the  average  line  Intensity  and  average  ^Ine  spacing*  ^y 
applying  the  random  Elsasser  oiodel  to  the  fund^utal  and  super¬ 
posed  higher  order  bands,  closed-form  solution^  are  obtained  for 
the  emlsslvlty  as  a  function  of  u  for  certain  :|imltlng  cases* 
Experimental  data,  vhere  available,  are  used  f(t>r  the  line  widths 
and  total  absolution  of  the  bands*  The  hermon:|lc  oscillator 
approximation  Is  used  to  estimate  strengths  of  higher  order  bands 
for  vhlch  experimented,  data  are  not  available*  This  analysis  Is 
applied  to  HC/,  HF,  CO,  and  NO  for  temperature^  ranging  from 

500  to  7000’K. 

I.  lUTRODCCTION 

A  first  approximation  to  the  Infrared  spect red.  emlsslvlty  of  a  diatomic 
gas  ceui  be  made  by  assuming  a  heumonlc  oscillator  mj>del  for  the  molecule*  * 
This  has  tho  advantage  of  simplicity,  but  the  very  pearly  synmetrlc  Inten¬ 
sity  distribution  vhlch  It  yields  does  not  closely  ^tch  the  asymnetrlc 
shape  vhlch  Is  found  experimentally  or  ced.culated  b;|r  a  llne-by-llne 
method*  ^ 


ICnSKSSUiii'UliM 
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11. 


Tc  is^rove  the  approximation,  the  BK>del  Is  assumed  to  be  an  anhar- 
Bonle  oscillator  vlth  the  first  approximation  to  the  vlbratlcn-rotatlon 
Interaction.  This  vlU  Improve  the  approximation  to  the  actxial  spectral 
emlsslvlty  without  resort  to  the  more  detailed  approach  of  considering 
the  emission  from  each  Individual  spectral  line. 


JEVELOPMEWT 

The  energy  levels  of  an  anhamonlc  oscillator.  In  the  first  approxi¬ 
mation  to  the  vibration-rotation  Interaction,  are  given  by 

3(V,J)  -  «e  ('  +5)  -  *,  ('  ♦  iT  +  y,  (»  +57  +  *.  V' 

♦Sg  J(J  +1)  -  0,  (»  +!)•  W 

(Refer  to  the  treatment  of  Stull  and  Plass^  for  notation  and  basic 
equations.)  The  frequency  of  a  transition  for  |iar|  «  1  Is  then  given  by 

M  *  E(v  +  1,  J  ')  -  E(v,J) 

-My  +B^(r(r  +1)  -  J(J  +1)1 

-  a,t(v  +|)  J'tr  +1)  -  (V  +|)  JU  +1)1,  (2) 

where 

Wy  =  Ug-2(v+l)  Ug  Xg+  (5(v+l)®+ J]  yg+ [  4(v+l)®+(v+l)  1  z^.  (3) 

In  tenns  of  ra(m  =  J  +  1  for  the  R-branch,  a  =  -J  for  the  P-branch), 

w  =  id  +  2B  a  -  a  [a (a  +  1)  +  2(v  +5)  a).  (4) 

This  equation  expresses  u  as  a  quadratic  function  of  a.  On  solving  for 
a,  this  becomes 


B 


a  = 


otgCv  +  1)  \  -  gg(v  +  l)f  -  ggCu  - 


My) 


(5) 
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(6) 


Note  that  Equation  (4)  iopllea  a  tsaxlcua  value  for  u  of 


B,  -  ge(v  *  1) 


Since  (B  •  a  (v  -»■  1)  ]/a  >  o  for  all  tsoleculee  under  consideration,  this 
Bodel  vlU  always  have  a  hand  head  in  the  R-branch.  Hence,  the  upper 
sign  before  the  radical  in  Equation  (5)  refers  to  the  aain  portion  of  the 
band  and  the  lower  sign  to  the  returning  R-branch.  (This  convention  for 
the  7  or  i  sign  vlU  be  used  throughout  this  analysis.) 

The  average  line  spacing  in  the  region  of  a  given  b  is 

a  ■  2|b^  -  Q^(v  +1  +  a) |.  (8) 

This  becoaes  [eiqpressed  in  terns  of  u  by  use  of  Equation  (!^)]» 

d(u)  «  2  -  o^(v  +  1))*  -  ci^(w  -  (9) 

for  the  entire  band,  Including  the  returning  R-branch. 

The  energy  tern  E^CJ)  (Equation  (6),  Stull  and  Plass^)  in  terns  of 
the  transition  frequency  «,  is  given  by 


EyO) 


-  {4®.  -  “e]  [®e 

[sg  -  (V  + 1)  -  ag(u)  -  U^)] 


Bg  -  (V  +  1) 


(V  +1)  0^ 


Cv 


Several  approxluntions  are  used  lu  the  nnalysis.  The  factor 
defined  by  Equation  (5)^  io  approximated  by 


”  =  “v  t  -  “P  (■  ^)]  • 
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The  Herm&n  and  Wallis^^  factor 

^  J 


la  approximated  by 


F  =  1  +  c 


-  ag(v  + 1)  T>/'  -  o^(v  +  l)f  -  a^(G)  -  My) 


] 


(12) 


vhere 


c 


.  8B^6 


(«) 


(9  Is  the  ratio  of  the  zero^^  to  the  first  order  term  In  the  expansion 
of  the  electric  dipole  manent).  The  partition  functions  Cp(T)  and  C(T) 
are  approximated  by 


Q  /<p\  -  ^ 

>  B  he 
e 


(Ilf) 


and 


G(T)  = 


kT 

B  he 
e 


exp  t.E(0.01  ho/KTl 
1  -  exp  (-hcm^kT) 


(15) 


Stull  and  Plass^  give  the  following  equation  for  the  integrated 
absorption  of  a  vibration-rotation  line: 


/  /  / 
V  V  J 


''  J  Gj(y,T)  0-J  I-  ''  J 

[•)  ®j-i, J '  *  (J  ^  'j+1, J ']  [■  -  ^  (a  i  s)]  • 


Since 


Vibj'"'  ^j+i,r " ^ 


8 


and  by  use  of  the  preceding  approxioationsi  this  becoaes 


kT 


u 


®e  *  l®e  f  -  0^(w 


‘X 

K 


EScperlmentol  data  are  used,  where  available,  for  the  total  band 
absorption  A  convenient  approxlantion  is  cade  to  Stull  and  Plaas*' 

Equation  (8);  naiaely, 


[|k(t)  -  E{0,0)J^  ^  -  0  ]  ' 

(18) 


SO  that  can  be  computed  for  any  teoperature  If  it  is  known  for  one 
tesperature. 

If  data  ore  not  available  for  higher-order  bands,  by  asBumlng  the 
harmonic  oscillator  approximation  and  using  Stull  and  Plass*'  Equations  (4) 
and  (9),  one  obtains 


• 

1 - 1 

(19) 

In  the  weak  line  approximation,  the  emlssivlty  Is  given  by 

r  ~\ 

€  =  1  -  exp  -  piS((A!)/d(u) 

0)  L  ' 

(20) 

S(c.)  =  Js^'-'-lc^) 

s 

(21) 

V  " 
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(The  SUB  over  v  Includes  all  hl^er-order  bands  which  contribute  signifi¬ 
cantly  to  the  eaission  In  the  region  of  the  fundamental,  including  bands 
associated  with  isotopic  species.) 

An  additional  approxlsation  Is  aade  in  th'  strong  line  approxlBation; 


In  general,  this  approxloatlon  causes  a  change  in  the  eBlssivity  of  less 

yf,,  (It  gives  the  sane  result  as  for  the  statistical  aodel. ) 

In  the  strong  line  approxiaatlon, 

-  1  -  exp  [.2a/  (p®/)*  B*  (u)/d(u)]  (22) 

where 

^  («)  -  £  [  («))*  +  (u))*]  ,  (23) 

and  a  is  the  average  half-width  of  the  line  at  l^tm  pressure. 

0 

A  detailed  discussion  of  the  range  of  validity  of  the  various  approxi¬ 
mations  is  given  by  Plass.*  For  exajsple,  Plass'®  Figure  2  indicates  that 
for  pure  CO,  when  *  0.6,  the  weak  line  approximation  is  accurate 
within  1036  when  0  >  5.  (This  requires  the  pressure  to  be  above  some 

valuc  whlch  Is  largest  (p  **  50  atm)  when  T  1200*K. )  The  strong 
line  approximation  is  accurate  within  lOjt  when  0  <  1.  (For  T  *  1200*K 
this  implies  a  maximum  p  1?  atm. )  A  lower  limit  to  the  validity  of  the 
strong  line  approximation  is  imposed  by  the  presence  of  Doppler  broadening, 
which  has  been  neglected.  The  Doppler  and  Lorentz  widths  are  conqjarable 
for  p  «  1/2  atm  at  T  «  1200*K. 

For  pure  HF,  the  strong  line  approximation  may  be  used  for  pressures 
up  to  a  certain  value  (a  maximum  p  «  55  atm  for  T  »  5000*K)  and  above  a 
minimum  determined  by  the  Doppler  broadening.  (The  Doppler  and  lorentz 
widths  are  comparable  for  p  »  atm  and  T  5000°K. ) 
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III.  DISCUSSION 


Results  of  this  analysis  arc  presented  (Figures  1  to  ^2)  for  SCI, 

SF,  CO,  and  NO  for  tesqperatures  ranging  fro&  300  to  7000*IC.  Ihe  presence 
of  isotopes  other  than  Cl^^  and  Cl^’’^  vas  Ignored.  Baission  froa  as  aany 
as  20  hands  vas  considered. 

Data  of  Benedict,  et  al.,^  for  HCi  were  used  to  sake  a  direct  coa- 
parison  with  the  nore  elaborate  calculations  of  Stull  and  Flass.^  Their 
value  of  130  i  7  ca*^  atm"^  for  the  fundamental  band  strength  is  probably 
too  low  (cf,  Babrov,  et  al.,*  I50  ±  5  cm“®  ato"^;  Weber  and  Fenner,® 

150  cm"®  atm"*’). 

Benedict  and  Flyler's®  value  of  260  cm"*  atm"*"  was  used  for  the  CO 
fundamental  band  strength  (cf.  Weber  and  Fenner,®  237  cm"*  atm"^). 

Kuipers'  data*^  were  used  for  HF,  and  Weber  and  Fenner's®  for  NO. 

Average  line  widths  were  cmiputed  as  described  by  Stull  and  Flass.  ^ 

It  should  be  noted  that  the  supplementary  canputation  made  by  these 
authors  for  emission  from  line  wings  in  the  frequency  range  beyond  a 
band  head  is  omitted  here. 

From  20  to  30  points  were  computed  for  each  curve.  Each  set  of  four 
curves  represents  a  charge  of  between  0.01  and  0.02  hour's  time  on  an 
IBM  704  digital  computer. 

A  conpairison  is  made  (Figures  35  34)  with  other  published  compu*' 

tatlons  of  diatomic  gas  emisslvltles.  Figure  33  ccagpares  Equation  (20) 
to  the  results  of  one  particular  computation  by  Stull  and  Flass^  (for  ECl 
at  2400"K  in  the  weak  line  approximation  for  pi  =  100  atm-cm)  in  which 
they  considered  the  emission  from  approximately  29OO  individual  spectral 
lines.  Fenner  and  Gray®  recently  published  a  comparative  conqnitatlon 
based  on  a  simple  harmonic  oscillator  model,  which  is  also  shown  in 
Figure  35 • 

The  curve  representing  Equation  (20)  follows  quite  closely  the  curve 
computed  by  Stull  and  Flass,  ^  which  is  presumably  the  most  accurate 
available  theoretical  computation.  The  agreement  would  be  slightly  closer 
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FIGURE  I  SPECTRAL  EMISSIVITY  OF  CO  AT  T  *  300*K 
(WEAK  LINE  APPROXIMATION) 
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FIGURE  3  SPECTRAL  EMISSIVITY  OF  CO  AT  T  •  I200*K 
(WEAK  LINE  APPROXIMATION) 
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FIGURE  4  SPECTRAL  EMSSIVITY  OF  CO  AT  T-IOOO*K 
(WEAK  LINE  APPROXIMATION) 


FIGURE  5  SPECTRAL  EMISSIVITY  OF  CO  AT  T>2400*K 
(WEAK  LINE  APPROXIMATION) 
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FIGURE  6  SPECTRAL  EMISSIVITY  OF  CO  AT  T»3000*K 
(WEAK  LINE  APPROXIMATION) 


FIGURE  7  SPECTRAL  EMIS51VITV  OF.  CO  AT  T*5000*K 
(WEAK  LINE  APPROXIMATION) 


FIGURE  9  SPECTRAL  EMISSIVITY  OF  NO  AT  T<300*K 
(WEAK  LINE  APPROXIMATION) 
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FIGURE  13  SPECTRAL  EMISSIVITY  OF  NO  AT  T°2400*K 
(WEAK  LINE  APPROXIMATION) 
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FIGURE  14  SPECTRAL  EMISSIVITY  OF  NO  AT  T<3O0O*K 
(WEAK  LINE  APPROXIMATION) 


FWUBE  <6  SPECTRAL  CliBSSIVITY  OF  NO  AT  T-TOOO*K 
(WEAK  LINE  APPROXIMATION) 


FIGURE  17  SPECTRAL  EMISSIVITY  OF  HCA  AT  T  >  300*K 
(STRONG  LINE  APPROXIMATION) 
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FIGURE  18  SPECTRAL  EMISSIVITY  OF  HCi  AT  T«600*K 
(STRONG  LINE  APPROXIMATION) 


FI6URE  21  SPECTRAL  EMISSIVITY  OF  HC^AT  T*  2400*K 
(STRONG  LINE  APPROXIMATION) 
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FIGURE  22  SPECTRAL  EWSSIVITY  OF  HC^AT  T* 
(STRONG  LINE  APPROXIMATION) 
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FI6URE  23  SPECTRAL  EMISSIVITY  OF  HC<  AT  T<5000*K 
(STRONG  LINE  APPROXIMATION) 
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FIGURE  27  SPECTRAL  EMISSIVITY  OF  HF  AT  T<I200*K 
(STRONG  LINE  APPROXIMATION) 


FIGURE  28  SPECTRAL  EMISSIVITY  OF  HF  AT  T«I800*K 
(STRONG  LINE  APPROXIMATION) 


FIGURE  29  SPECTRAL  EMISSIVITV  OF  HF  AT  T*2400*K 
(STRONG  LINE  APPROXIMATION) 


FIGURE  30  SPECTRAL  EMISSIVITY  OF  HF  AT  T>3000*K 
(STRONG  LINE  APPROXIMATION) 


FIGURE  32  SPECTRAL  EMISSIVITY  OF  HF  AT  T«7000*K 
(STRONG  LINE  APPROXIMATION) 


FIGURE  33  COMPARISON  WITH  OTHER  PUBLISHED  DATA  OF  SPECTRAL  EMISSIVITY 
COMPUTED  FOR  HCIAT  T«2400*K  IN  THE  WEAK  LINE  APPROXIMATION 
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FIGURE  34  COMMRISON  WITH  OTHER  PUBLISHED  DATA  OF  SPECTRAL  EMISSIVITY 
COMPUTED  FOR  NO  AT  T>3000‘K  M  THE  WEAK  LINE  APPROXIMATION 


if  tl..'  previously  mentioned  correction  for  line  wing  emission  beyond  the 
bond  head  had  been  it  eluded  in  this  analysis.  Pi’esuraably  the  tail  drawn 
by  Stull  and  Plass^  in  the  i*egion  near  5500  cm"^  was  computed  for  one 
particulai*  combination  of  path  length  and  pressure,  pertiaps  for  p  =  100  ata 

and  '2=1  cm. 

Figure  34  shows  the  results  of  a  computation  by  Breene^°  (for  HO  in 
air  at  3000‘*K,  density  ratio  =  10,  thickness  =  5  cm)*  Penner,  Sulzmann, 
and  Ludwig's®  r^'^ults  for  the  simple  harmonic  oscillator  model  are  also 
shown.  It  sho^tLa  be  noted  that  for  comparison  all  three  of  these  compu¬ 
tations  utlli-sed  Brecne'6^°  multiplicative  factor  of  1.88  for  the  average 
path  length  (l.e.,  for  a  5-cm-thlck  layer,  a  value  of  f  =  9.4  cm  was  used). 

A  glance  at  Figure  34  shows  that  Equation  (20)  is  not  much  more  suc- 
cesEi'ul  than  the  simple  harmonic  oscillator  approximation  In  matching 
Breene's^°  computation.  The  correction  for  line  wing  emission  would  add 
a  tall  to  the  curve  ripresentlng  Equation  (20)  and  extend  it  somewhat 
farther  past  2000  cm“^.  (Since  the  line  widths  are  of  the  order  of  2  cm"^, 
this  addition  would  not  be  very  significant.)  The  highly  aBynnetrlc  peaks 
in  emisslvlty  indicated  by  Equation  (20)  seem  more  reasonable  than  the 
nearly  symmetric  peaks  shown  by  Breene.  (Note,  in  conqparlson,  that  for 
HC2,  the  inclusion  of  the  Herman  and  Wallls^^  F-factor  has  caused  a 
considerable  augmentation  of  emisslvlty  at  lower  frequencies  and  diminution 
at  higher  frequencies  which  makes  the  peaks  nearly  equal.  This  factor  is 
negligible  for  the  more  symmetric  NO  molecule  so  that  considerably  more 
asymmetry  would  be  expected  on  this  basis. ) 

Emlsslvltles  have  been  computed  (but  not  presented  here)  for  HCi  and 
HF  in  the  weak  line  approximation  and  for  CO  and  NO  in  the  strong  line 
approximation  for  the  temperatures  considered  here.  Also  for  several 
intermediate  temperatures,  emlssivlties  have  been  computed  for  all  four 
molecules  in  both  approximations.  •  It  is  planned  to  extend  the  computa¬ 
tions  to  other  molecules  and  also  to  consider  Doppler  broadening  as  well 

as  Lorentz  broadenin,:. 
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